arXiv:astro-ph/0503046vl 2 Mar 2005 


Mon. Not. R. Astron. Soc. 000. [111311 0 Printed 2 February 2008 (MN IATJ 3 X style file v2.2) 


Simulation of a population of isolated neutron stars 
evolving through the emission of gravitational waves 

C. Palomba 1 * 

1 Istituto Nazionale di Fisica Nucleare, sezione di Roma 
P.le A. Moro, 2 - 00185 Roma, Raly 


ABSTRACT 

We study, via a Monte Carlo simulation, a population of isolated asymmetric 
neutron stars where the magnitude of the magnetic field is low enough so 
that the dynamical evolution is dominated by the emission of gravitational 
waves. A starting population, with age uniformly distributed back to 100 
Myr (or 500 Myr) and endowed with a birth kick velocity, is evolved in the 
Galactic gravitational potential to the present time. In describing the initial 
spatial distribution, the Gould Belt, with an enhanced neutron star formation 
rate, is taken into account. Different models for the initial period distribution 
are considered. The star ellipticity, measuring the amount of deformation, is 
drawn from an exponential distribution. We estimate the detectability of the 
emitted gravitational signals by the first and planned second generation of in¬ 
terferometric detectors. Results are parametrized by the fraction of the whole 
galactic neutron star population made of this kind of sources. Some possible 
mechanisms, which would make possible the existence of such a population, 
are discussed. A comparison of the gravitational spin-down with the braking 
due to a possible interaction of the neutron star with the interstellar medium 
is also presented. 
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1 INTRODUCTION 


Continuous gravitational signals emitted by rotating asymmetric neutro n stars are p romising 


wave de 


ectors like GEO (Willkeetal 


2002|) LIGO 


sources for interferometric gravitational 

( Ambramovici et ahlll992 h TAMA300 ( Tsubonolll995tl and VIRGO (jCaron et al .11 19971) . The 


detectability of a source, emitting due to a deviation of its shape from axi-symmetry, basically 
depends on the detector sensitivity, on the source rotation frequency, distance and ellipticity. 
Not less important for the feasibility of the search are, however, also the source position in 
the sky and spin-down rate. Up to now about 1500 neutron stars, most of which are pulsars, 
have been detected in the Galaxy through their electromagnetic emission. Standard theory 
of stellar evolution suggests that about 10 9 neutron stars should exist in the Galaxy and, 
of these, about 10 5 should be active pulsars. We detect only a small fraction of this subset 
because either the radio beam does not intersect our visual line or the intensity of the 
emission is too low. In fact, the sample of detected pulsars is increasing as the sensitivity of 
the surveys increases. For most of these stars the rotation frequency, the spin-down rate, the 
location in the sky and the distance are known, while the ellipticity is completely uncertain. 
This makes difficult to get reliable estimates of their gravitational emission. An upper limit 
to the ellipticity, and then to the amplitude of the emitted gravitational signals, is tipically 
given assuming that all the observed spin-down is due to the emission of gravitational waves. 


In this case e ^ 1.8 • 10 5 V PP 3 , being P , P respectively the rotation period and its first time 
derivative and, for most normal pulsars, values of the order of 10 -3 -ylCV 2 are obtained. Using 
this assumption, in Fig. 0 we show the effective amplitude of the gravitational wave signals 
emitted by the presently known isolated pulsars 1 , after using optimal filtering techniques 
and four months of observation time, versus the target sensitivit y of the Virgo detector 


and a possible sensitivity curve of an advanced Virgo detector (Punturo 


2004). Similar 


plots can be obtained for the LIGO interferometers. It comes out that very few observed 
pulsars would be on the verge of detectability, for the first generation of interferometers. 
Tens of pulsars would be detected by an advanced detector. On the other hand, as it is 
well known, it is quite likely that true cllipticities are much lower than the upper limits. 


According to 


Ushomirskv. Cutler & Bildsten (2000) the maximum ellipticity is probably 


1 http://www.atnf.csiro.au/research/pulsar/psrcat/. We do not take into account pulsars in binary systems for which the 
detection would be complicated by the system orbital Doppler effect and by the stochastic fluctuations of the rotation period 
due to the accretion of matter. 
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Figure 1. Effective amplitude of the gravitational wave signals emitted by known isolated pulsars versus the target sensitivity 
of the Virgo detector (continuous line) and a possible sensitivity curve for an advanced Virgo detector (dashed line). The total 
observation time is four months and the use of an optimal data analysis method is assumed. The observed spin-down is assumed 
to be due completely to the emission of gravitational radiation. 



Figure 2. Same as Fig.(l) but taking the source ellipticity as the minimum between 10 6 and that obtained assuming the 
observed spin-down is completely due to the emission of gravitational radiation. 


< 10 6 . This makes rather unlikely the detection of gravitational waves from the observed 
pulsar population: Fig. © is the same as Fig. m but taking the ellipticity as the minimum 
between 10~ 6 and the value obtained assuming only gravitational spin-down. The observed 
spin-down of millisecond pulsars implies a maximum ellipticity less than 1CB 8 , which is a 
plausible value. In this case we see that no detection is expected for the first generation of 
interferometers while few sources (4 at present), all millisecond recycled pulsars, could be 
detected by an advanced detector. By the way, also this is just an upper limit. Nevertheless, 
the search for gravitational signals from known neutron stars is an important task because, 
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even if no detection should come, we can at least put upper limits on the neutron stars 
ellipticity, thus increasing our knowledge on the neutron star structure. 

An alternative possibility is the existence of a population of asymmetric rotating neutron 
stars with magnetic held low enough so that the spin-down is dominated by the emission of 
gravitational waves. For such objects, as the electromagnetic emission does not contribute 
much or at all to the spin-down, a larger amount of gravitational radiation would be emitted 
in a given period of observation, in comparison to a standard pulsar with the same degree of 
deformation. Obviously, it would be very difficult to detect this kind of neutron stars in the 
electromagnetic band, unless very near and young (through their thermal emission) or very 
near and interacting with the interstellar medium. Hence the necessity of performing a blind 
search in the data of gravitational detectors, i.e. without strong assumptions on the position 
of the objects, on their rotation frequency and spin-down rate. The existence of such class 
of stars, which for simplicity we call GWDNS (Gravitational Wave Driven Neutron Stars) 
or gravitars , cannot be excluded both on observational and theoretical basis. For instance, 
the lack of a one-to-one correspondence between observed supernovae and observed compact 
remnants leaves room for a non negligible formation rate for GWDNS. In this context, the 
presence of a low magnetic held neutron star has been recently suggested for the remnants of 
SN 1987 A 


GWDNS could be a subset of pulsars born with a very 
low magnetic held or the product of a different evolutionary channel. Different mechanisms, 
like magnetic axis alignement, the freezing of the magnetic held in a newborn neutron star 
or magnetic held decay, could permit the existence of GWDNS. 

In this work we perform a Monte Carlo simulation of a possible GWDNS population, 
estimating the number of expected detections with the hrst and planned second gener¬ 
ation of gravitational interferometers. Several simulations of the neutron star or pulsar 


(Paczvnski 

1990*; 

Na.ravan & Ostri 

cer 

1990; 

Lvne & Lor 

mer 

1994; Hansen & Phinnev 

1997 

Arzoumanian. Chcrnoff & CordeJ 

2002; 

Gonthier et al. 

2002 

b Only in few cases the detec- 

tion of gravitational waves has been discussed ( 

e.Ef. Reeimbau & de Freitas Pacheco 

2003) 


The existence of a population of neutron stars whose dynamical evo 


the emission of gravitational waves has been addressed by 


Blandfordl (jl987f) but, as far as we 


ution is dominated by 


know, no attempt has been done, up to now, to make a quantitative study of the problem. 
We do not have an observed population against which the simulated one can be checked 
so that different choices for the parameters are done and the dependency of the results on 
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them is discussed. The initial position and kick velocity are described by the probability 
distributions typically used to simulate the standard pulsar population. The contribution to 
the total neutron star formation rate due to the Gould Belt is taken into account and three 
different distributions for the initial rotation period are considered. The mean value of the 
cllipticity, distributed exponentially, is taken between 1CT 8 and 1CT 6 and the abundance of 
GWDNS, expressed as fraction of the whole galactic neutron star population, is considered 
as a free parameter. Given the initial spatial, initial period and birth kick velocity distri¬ 
butions, GWDNS are evolved solving the equations of motion in the galactic gravitational 
potential, assumed axisymmetric for simplicity. The stars age in drawn from a uniform dis¬ 
tribution back to 100 Myr (or 500 Myr, when advanced detectors are considered). Once 
the population has been evolved to the present, the detectability of the emitted continuous 
gravitational signals is estimated, for the first and planned second generation of interfero¬ 
metric detectors, assuming that a blind search is done. This is mandatory because, contrary 
to what happens in targeted search where the source parameters are well constrained, here 
we need to explore a large portion of the source parametrer space. This implies that optimal 
data analysis methods cannot be applied, due to the unrealistic amount of computing re¬ 
sources they would require. Non optimal procedures, with reduced sensitivity but allowing 
to strongly cut the computational weigth of the analysis, are therefore needed. 

The plan of the paper is the following. In Section 2 we review the main relations con¬ 
cerning the time evolution of neutron stars rotation period due to the presence of a dipolar 
magnetic held or of an asymmetry in the star shape. In Section 3 the sensitivity of a blind 
search for the gravitational waves emitted by asymmetric rotating neutron stars is given. In 
Section 4 the Monte Carlo simulation model is presented and in Section 5 the results of the 
simulations are discussed. In Section 6 we describe some possible mechanisms which could 
leave space to a non negligible fraction of GWDNS. In Section 7 conclusions are drawn. 
Finally, in the Appendix we compare the gravitational spin-down of a neutron star with the 
spin-down due to its possible interaction with the interstellar medium. 


2 SPIN EVOLUTION BASIC RELATIONS 

The spin-down rate of a neutron star can be expressed, in the general case, as 
H = - £ K,Q rH (1) 

i 
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where f2 is the rotation pulsation, n* are the braking indexes which depend on the spin-down 
mechanisms at work and Ki are constants depending on the neutron star structure and 
characteristics. In the case of simple electromagnetic dipolar emission Eq. m becomes 


2 B^ 

^dip J ^ ^ ^ 


( 2 ) 


where B is the magnitude of the dipolar magnetic held, a is the angle between the magnetic 
axis and the rotation axis, R is the neutron star radius and / is the star momentum of inertia. 
If the evolution is dominated by the emission of gravitational waves in a non-axysimmetric 
neutron star, we have 


h gw = ~^^lle 2 n 5 


5 c 5 


( 3 ) 


where e = 11 1 * 2 is the star cllipticity, being Ii.R.G the star momenta of inertia along 


_ h-h 
h 

three principal axes. From this equation the rotation frequency as a function of time can be 
expressed as 

2048tt 4 G 


m = k + 


rIGt 


( 4 ) 


5 c J 

where Po is the initial rotation period. The final frequency is independent on the initial 
period if the inequality 


P 0 < 0.0443 


I 


1.12 • 10 38 kg m 2 J V10 ~ 6 J V 1 G V r 


t 


is satisfied and, in such a case, the rotation frequency at present is 


m - 23 


I 


1.12 • 10 38 kg m 2 


10 


-6 


t 


Hz 


( 5 ) 


( 6 ) 


\1 Gyr J 

We note that even very old and rather highly deformed GWDNS can continue to spin to 
a rather high frequency at present. GWDNS only slightly deformed can maintain nearly 
unchanged their rotation frequency along their whole life. 

We are interested in comparing the contributions to the spin-down rate due to the emis¬ 
sion of gravitational waves and to the presence of a dipolar magnetic held. In this way we 


can establish whic 


i_condition the magnetic held must satisfy in order to have a GWDNS. 


Let us introduce (Palomba 


2000) the ratio among the gravitational and the dipolar held 
induced spin-down rates, expressed by Eqs. (EED: 


Y = 


n 


gw 


n 


= 2.3-10 


-7 


I 


dip 


1.12-10 38 kg m 2 I V 10~ 6 


2 / f / B \“ 2 / R W 6 


/ 


100 Hz, 


\10 12 GV V10 km) 


sm a 


-2 


The condition Y > 1 corresponds to 
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B < 4.75 • HE 


I 


^c) ( ~ ) (~i—) 3 sin a” 1 G 

10~ 6 / l 100 Hz \10kmJ 


( 8 ) 


1.12 ■ 10 38 kg m 2 

Obviously, this relation is less stringent for high values of the neutron star cllipticity and 
frequency and small values of the angle a. In order to have a GWDNS is not sufficient that 
the condition Y > 1 is satisfied at the birth. Due to the spin-down, the ratio Y, even if 
> 1 at the beginning, decreases in time and could go below 1, unless the magnetic held 
decays or the magnetic axis aligns to the rotation axis. These possibilities will be discussed 
in Section 6. Assuming that the ”non-gravitational” spin-down is described by Eg. ( 121 ). we are 


implicitly considering a neutron star in the so-called ejector phase (|Popov et al 


2000). This 


is the evolutionary stage to which pulsars belong. Note, however, that not all ejectors are 
pulsars. Two other phases are foreseen for neutron stars having small values of the magnetic 
held and/or long rotation periods: propeller and accretor. In these phases an interaction of 
the star with the sorrounding interstellar medium takes place and in the Appendix we will 
see that, unless the magnetic held is nearly vanishing, the spin-down rate in these stages is 
dominated, for e > 1CD 9 , by the emission of gravitational waves. 


3 DETECTION OF CONTINUOUS GRAVITATIONAL SIGNALS 

The gravitational wave signal emitted by an asymmetric neutron star, rotating around one of 
its principal axis of inertia, has a frequency twice the rotation frequency and the amplitude 
at the detector, assuming optimal orientation, can be expressed as 


hn = 1.05 • 10” 27 


1.12 • 10 38 kg m 2 J V10 ~ 6 ) \ 100 Hz 


6 \ f fgw 
10-V V 100 H 


10 kpc s 


(9) 


where f gw — is the gravitational wave frequency and r the distance from the detector. 
From Eq. ©, assuming e ^ 10 6 , most GWDNS emit gravitational waves with frequency 
> 10 Hz, independently of their age, i.e. within the sensitivity band of interferometric 
detectors. However, we will see in Section 0 that almost all GWDNS producing signals 
detectable by the hrst generation of interferometric detectors are younger than 10 Myr. 
The expected amplitude of the gravitational signals is so low that they will be typically 
buried in the detector noise. Then, very long integration times (of the order of months) 
and suitable data analysis techniques are needed to extract them. It is important to take 
into account that the signal ’’seen” by the detector is not really monochromatic due to 
the source spin-down and the Doppler effect associated to the detector motion. Moreover, 
the search for GWDNS must be necessarily blind , that is it must be performed on a large 
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portion of the source parameter space, because we cannot put strong constraints on their 
position in the sky, on their rotation frequency and on their spin-down rate. This means, 
in practice, that we want to search sources everywhere on the celestial sphere ( full-sky 
search), with signal frequencies from the low frequency cut-off of interferometers up to 
about 2 kHz , and a spin-down rate as large as possible, compatibly with the available 
computing resources. Such kind of analysis cannot be performed with the classical optimal 
methods used in gravitational wav e experiments, based on the so-called match ed filter, which 


can be used for targeted searches (je.g. Astone et al 


2QU2; 


Abbot et al 


2004J), but that are 


unfeasible for wide area searches, due to the enormous computing power they would require. 


Alternative, non op 

:ima 

methods have been developed by various groups in the world 

(Frasca & La Posta 

1991; 

Astone. Frasca & Papa 

1997; 

Bradv et al. 

1998; 

Schutz & Papa 


199fll 


Frasca 


2000). They reach a lower level of sensitivity but strongly cut the needed 


computational power. All these methods are based on some kin<Aof alternation of coherent 


and incoherent steps, t hese l ast ty pical! 


transform (Frjiggg 


our group (iFrasca 


2000 


2000 ; 


[to 


Krishnan et 


Brocco et al 


al 


jased on the Radon ( Bradv et al 


1998 ) or Hough 


2004). In the h ierarchical method developed by 


2003 


Frasca & Palomba 


12004 s ) 


smallest detectable signal, with signal-to-noise ratio ~ 2, is 

4 sl(f) 


hmin(f) 


P 


the amplitude of the 


( 10 ) 


TobsTFFT J 

where Sh(f) is the detector noise power density, measured in Hz~ l , T obs is the total ob¬ 
servation time, and T FF t is the time duration of the short FFTs computed from the data 
at the beginning. T F ft varies between 1048 s and 16667 s, depending on the frequency 
band which is analyzed. The loss factor (3 is a number around 2, which exact value depends 
on the analyzed frequency band, on the number of candidates which are selected after the 
first incoherent step and on the assumed minimum spin-down decay time. The larger is the 
available computing power and the wider can be the search. A computing power of about 
1 Teraflops will allow to perform a full-sky analysis of 4 months of data, searching sources 
emitting signals with frequencies up to 2 kHz and with a spin-down time r = j > 10 4 yr. 

From Eqs. dump we can derive the maximum distance r max at which a source, emitting 
gravitational waves at a given frequency and with a given cllipticity, can be detected. In 
Fig.© r max is plotted as a function of the frequency for different values of the cllipticity. 
A total observation time of 4 months is considered and the target Virgo sensitivity is used. 
The discontinuities are due to the fact that different values of the FFT duration, Tfft, 
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Figure 3. Maximum distance at which the gravitational signal emitted by an asymmetric rotating neutron star could be 
detected, as a function of the signal frequency, assuming the hierarchical procedure developed by our group is applied and 
using the Virgo target sensitivity. The search parameters are: total observation time of four months, minimum spin-down decay 
time of 10 4 yr , 10 9 candidates selected after the first incoherent step. The different curves corresponds to different values of 
the star ellipticity, starting from 10 —6 (upper curve) to 10 —9 (lower curve). 


are used for different frequency bands. Looking at the figure we see that if e ^ 10 6 and if 
the maximum rotation frequency we search for is 1 kHz , the maximum distance at which a 
source can be detected, in a blind search, is about 6 kpc. On the other side, assuming that no 
nentron star is located at a distance from the Earth less than 10 pc, the minimum ellipticity 
needed to have a detectable source (emitting at the high end of the allowed frequency range) 
is ~ 1.4 • 10“ 9 : no source can be detected, not even in principle, if its ellipticity is less than 
that value. If the typical ellipticity is in the range 10 “ 8 -r 10“' we would need to have sources 
located in the solar neighbourhood and rotating at high rate. Moreover, sources with rotation 
frequency less than ~ 25 Hz cannot be detected, indipendently of the ellipticity. In Section 
5 we will see that detected sources in the simulation have parameters in a narrower range. 


4 MONTE CARLO SIMULATION MODEL 

In this section we describe the model used to simulate the population of GWDNS in the 
Galaxy. As we do not have any observational clue on what the parameters distribution 
of GWDNS could be, we use for initial position and kick velocity the same distributions 
typically used for pulsars. Three different models are considered for the distribution of 
initial periods. The ellipticity follows an exponential distribution with mean value treated 
as a parameter. Obviously we cannot exclude that GWDNS have parameters distributed 
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in a completely different way. We take the distance of the Sun from the galactic centre as 
Rq = 8.5 kpc. 


4.1 Space distribution 

We ass um e that the initial space distribution is described by two exponentials in the radial 
distance from the galatic centre, R = \Jx 2 + y 1 , and in the height over the galactic plane z: 

p(R) = (11) 

rto 

with scale factor R 0 = 3.2 kpc, and 

p(z) = e~ y (12) 

with Zo = 0.075 kpc. The results of the simulations are rather insensitive to the choice of the 
scale factors. The initial right ascension, </>o, is chosen from a uniform distribution between 
0 and 27r. The angle is measured respect to the direction between the galactic centre and 
the Sun. With these choices we find that about 1100 neutron stars were born in the galactic 
disk up to a distance of 3 kpc from the Sun in the last 4 M yr (as suming a neutron star 
formation rate of 0.02 yr~ x ), in good agreement with Popov (|2Q04J ). The presence of the 
Gould Belt is also taken into account. It is described as a thin disk with radius of 300 pc and 

le galactic plane and centered at 100 pc from the 


thickness 60 pc, tilted of 18° respect to t 
Sun in the Galactic anticenter direction (Popov 2004). The Gould Belt is characterized by 


an overabundance of massive stars and its age is about 40 Myr. Following Popov ( 200411 . we 
assume a constant star formation rate inside the belt of 27 neutron star per Myr. To mimic 
the observed lack of massive stars in the central region of the Gould Belt, we assume that 
no neutron s tar wa s born in its central region up to a distance of 150 pc from the center 
Jpopov et alibofliiL We find that a fraction ~ 1.5 • 10~ 3 of the neutron star population with 
age less than 40 Myr was formed in the solar neighbourho od, at a distance from the Sun 


less than 600 pc, and this is in agreement with the result of 


Popov et al. 


(2003) 


4.2 Kick velocity distribution 

We assume that the distribution of the modulus of the kick vcloci 


: v Vhirh imparted to neutron 


stars at birth is that given by Arzoumanian. Chernoff & Cordes (2002): 


/ 


p(v k 


ick, 


= 4?rV lick 


Wi 


{ 27T <) 


-e 2<T «i + 


(i-wi) -'a* 


\ 


( 2yr <) 


v 2 


(13) 
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where w\ = 0.4, a Vl = 90 kms _1 , cr V2 = 500 kms~ l . The angles <p v ■ 0 V , defining the direction 
of the kick velocity vector, are drawn from uniform distributions between, respectively, [0, 27r) 
and [—f, f]- 


4.3 Initial period distribution 


We use three different models for the initial period distribution. The first one is a standard 
log-normal distribution 
1 


-(log Pp-log Pq ) 2 


p(Po) = 


V / 27T <J\ogP Q Pq 


log Pq 


with logP 0 = —2.3 and ai og p 0 = 0.3 (Arzoumanian, Chernoff & Corde; 


JH). 


(14) 


GWDNS with 


Po < 0.5 ms are excluded from the simulation. The second model is like the first one but we 


set to 10 ms all initial periods shorter than t 

lis value. In this way we try to mimic the possi- 

blc presence of excited r-modes ( 

Owen et al. 

1998 

) which could rapidly increase the rotation 


period of a fast rotating newborn neutron star. According to lAndersson. Kokkotas & Schutz 


), see section 3 of their paper, the final period of a newborn neutron star, once the 
r-modes have become stable, is between ~ 5 ms and ~ 15 ms, depending on the cooling 
scenario. We take 10 ms as a intermediate value. The third model consists of a uniform 
distribution of periods between 2 ms and 15 ms, mimicing in this way the combined effect 
of r-modes , which tend to increase the rotation period, and of the matter fall-back after 
supernova exp losion which, on the other hand, increases the star angular momentum. Ac¬ 
cording to Watts & Andersson (2002), the resulting period after one year or so, is largely 
independent on the period at birth and depends mainly on the neutron star magnetic field, 
with lower magnetic field and lower r-modes saturation amplitude resulting in lower periods 
(see section 4 of their paper and, in particular, figures 2 and 4). They find that neutron 
stars with magnetic field B ~ 10 12 G have final rotation periods, about one year after the 
supernova explosion, somewhere between ~ 2 ms and ~ 15 ms. For GWDNS we could ex¬ 
pect substantially lower values, but given all the uncertainties we prefer to take the period 
uniformly distributed in that range. 


4.4 Age distribution 

The age of GWDNS is chosen from a uniform distribution between 0 and 100 Myr or 
500 Myr. The reason for these choices is the following. We have computed the present fre¬ 
quency distribution of the signals emitted by a population of GWDNS, evolving according to 




















12 C. Palomba 



f [Hz] f [Hz] 

gw 1 1 gw 1 ‘ 


Figure 4. Final emission frequency distribution for a population of GWDNS with ages up to 10 Gyr , for different fixed values 
of the star ellipticity. 


Eq.(JH), for different values of the star ellipticity and assuming a maximum age t max = 10 Gyr. 
From Fig. 0 we see that most of the sources emit in the detector frequency bandwidth, even 
for high values of the ellipticity. Then, in principle, we should simulate the population of 
GWDNS with ages up to the age of the Galaxy. In practice, however, we have verified that all 
detected sources in the simulation have ages much lower than that. In the case of the initial 
Virgo detector, in particular, most detectable sources have age less than about 1 — 2 Myr. 
This is a consequence of the fact that mainly fast rotating (and then young) objects can pro¬ 
duce detectable signals. As already said in Section 3, detectable GWDNS typically belong 
to the solar neighbourhood (see Section 5.1 for a more quantitative discussion). A typical 
neutron star with a kick velocity of about 400 kms -1 will move, in that time, of about 
600 pc respect to the Sun. Then, considering also the over production of neutron stars in 
the Gould Belt, sources today in the solar neighbourhood are likely to be born in the solar 
neighbourhood. To be conservative, we have chosen t max = 100 Myr for the first generation 
of interferometric detectors and t max = 500 Myr for advanced detectors. Assuming a con¬ 
stant neutron star gslsctic formation rate of ~ 0.02 yr" 1 we have that N tot ~ 2 • 10 6 neutron 
stars have been formed in the last 100 Myr (or 10' in 500 Myr). The number of GWDNS 
will be a fraction A of N tot that we treat as a free parameter. N tot is the number of neutron 
stars we have used in the simulations. 
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4.5 Ellipticity distribution 

The typical ellipticity of neutron stars is largely unknown and also the allowed maximum 
value, depending on the properties of matter at very high densities, i s very un cert ain. At 
prese nt, maximum values e max ~ lCT'-hlCT 6 are considered plausible (Ushomirskv, Cutler & Bildsten 
201)0). We cope with our poor knowledge by using, for the ellipticity, the probability dis¬ 
tribution which satisfies the principle of maximum entropy , that is that leaves us with the 
largest uncertainty, without introducing any additional assumption or bias. The probability 
distribution p(x) maximizing the differential entropy H(x) = — f p(x) log p(x)dx, under the 

constraint that the expectation value is given, is the exponential distribution 
_ 

e T 


p(‘) = 


(15) 


_ emax \ 

1 — e r J 

where e max is the maximum possible value for e. The relation between r and the mean value 


of the distribution, e, is 

_ €max 

e = r 


e max 

e r — 1 


(16) 


In the simulation, chosen a value for the mean ellipticity e, the parameter r is calculated 
numerically using the previous equation. With the choice e max = 2.5 • ICC 6 , the maximum 
mean value is ICC 6 . In Section 5, we will report results, concerning the detection rate, for 
the different initial period distributions and for mean ellipticities e = ICC 6 , ICC 7 , ICC 8 . 


4.6 Galactic gravitational potential energy 

We calculate the motion of GWDNS in the Gala xy usi ng the following expressions for the 


galactic gravitational potential energy (jjPaczyhski 

\ GMi 


1990): 


\j R 2 + [oj + \J~[z 2 + 


(17) 


where i = 1,2 correspond respectively to the spheroid and the disk, and 
GM r / 1, r 2 1 - 2 

| 1 -I-In f 1 4- - 


®h(r) = 


I 1 ry*** ry* ny* 

1 + - In (1 + -|)-In (1 + —)-- arctan — 


rt 


r 


r. 


(18) 


which describes the galactic halo, with r = \Jx 2 + y 2 + z 2 . Following 


Gonthier et al. 


( 2002 ), 


we have added the constant term in the halo potential, which does not affect star trajectories 
but is important if one wants to evaluate the number of objects which escape from the 
Galaxy. The parameters are 

Mi = 1.12-10 10 M 0 , a\ = 0, bi = 0.277 kpc (19) 
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M 2 = 8.07 ■ 10 10 M q , a 2 = 3.7, b x = 0.20 kpc (20) 

M c = 5.0 ■ 10 10 Mq, r c = 6.0 kpc, r h = 41 kpc (21) 


From the galactic potential we calculate the rotation velocity as a function of the position 
(R, z ) (being the potential axisymmetric, the rotational curve does not depend on </>) using 
the relation 
d<f> 


Vrot R ()r 

with <h(i?, z) = (R, z) + <h 2 (R, z) + &h.(R, z). We hnd 

R 2 


v, 


rot, = GM , 


GM, 


R 2 + [at + \Jz 2 + b if] 2 

i? 2 / 


^rot,h 


R 2 + z 2 


1 - 


V^ 2 + 2 s 


arctan 


'VR 2 + z 2 ' 


( 22 ) 


(23) 


(24) 


The modulus of the total rotational velocity is v rot = \Jv 2 ot x + n 2 ot 2 + ^rof,/r O nce we have 
extracted a triple (vkick, 4>v> @v) for the kick velocity and a triple ( R, z, fo) for the star initial 
position, we integrate the equations of motion: 


R = 
with 


< 9 $ 


dR 


eff 

— ■ Z = 




eff 


dz 


(25) 


®eff(R, Z ) = $(R, Z) + 


2 R 2 


being L z the z-component of the star angular momentum which is a constant of motion with 
value depending on R, v rot , Vki C k, 4>v > do> ■ This system of two equations of the second order 

can be reduced to a system of four equations of the first order which we numerically solve 
with initial conditions R(t = 0), z(t = 0 ),u(t = 0) = R(t = 0), v(t = 0) = z(t = 0) using a 
hfth order Runge-Kutta method with adaptive step size. In the trajectory of each star the 
total energy is conserved to better than one part in 10 6 . 


5 RESULTS OF THE SIMULATIONS 

We have done several simulations of a population of N tot = 2 • 10 6 GWDNS using the 
probability distributions previously described. N tot is the number of neutron stars we expect 
to be born in the Galaxy in the last 100 Myr, assuming a constant formation rate of 
0.02 yr -1 , see the discussion in Section 4.4. In Table [T] the properties of the different models 
are summarized. In each run the number of detected sources is computed, assuming the 
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Table 1. Summary of model parameters. Concerning the initial period distribution function, p(Po), ’’standard” refers to the 
period distribution described by Eq. ITU : ”r -modes” corresponds to the same period distribution of Model 1 but in which all 
periods less than 10 ms are set to 10 ms, while ’’r-modes+fall-back” describes a uniform distribution between 2 ms and 15 ms, 
see Section EJfor a detailed description. 


Model 

p(Po) 

e 

1 

standard 

icr 6 

2 

standard 

io - 7 

3 

standard 

10“ 8 

4 

r-modes 

IO" 6 

5 

r-modes 

10- 7 

6 

r-modes 

icr 8 

7 

r-modes+fall-back 

icr 6 

8 

r-modes+fall-back 

IO" 7 

9 

r-modes+fall-back 

icr 8 


following search parameters, already discussed in Section 01 total observation time T+ = 4 
months, full-sky search, signal frequency up to 2 kHz , minimum spin-down time r = 10 4 yr, 
10 9 candidates selected after the first incoherent step. Some quantities related to the detected 
sourcers are also computed: the 50 th and 90 th percentiles of the signal frequency, distance, 
declination angle, and age distribution functions. To reduce statistical fluctuations results, 
for each model, are averaged among 300 independent runs for initial Virgo and over 50 
independent runs for the advanced Virgo detector. Each run takes about 2 hours on a 
Xeon 2.4 processor machine. We have performed the whole set of simulations within the 
Production Grid 2 of the italian Istituto Nazionale di Fisica Nucleare, consisting of a large 
set of geographically distributed computational and storage resources covering more than 
20 sites. This has largely reduced the time needed for the computations. 

Independently of the model, we have found that ~ 45% of the neutron stars population 
will escapejfom tiny Galaxy, due_to the high kick velocity. This is a well-known result found, 
e.g., by Arzouma nian. Chern o ff & C ordes (2002). 


According to 


Popov et al 


12003^ . the number of neutron stars within ~ 1 kpc from the 


Sun and with age less than ~ 4 Myr should be, at present, around 100. The number we 
find is ~ 140, comparable with that. Of these, ~ | were born in the Gould Belt. Most of 
the detected sources, with the first generation of interferometers, have ages less than that 
value. Then, the presence of the Gould Belt impacts in a significant way on the expected 
number of detections. 

A negligible fraction of the detected sources has a spin-down decay time less than 10 4 yr. 

Let us now discuss in more detail the results of the simulations considering separately 


2 http://grid-it.cnaf.infn.it/ 
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Table 2. Results for Virgo detector. The first column indicates the model, according to the classification of 
Tab, ffl : the second column shows the number of expected detections in units of A, which is the fraction of 
the total neutron star population made of GWDNS; f 50 ,^ 50 , \0\ 50,^50 and / 9 o,^ 90 , |0| 90,^90 give the values 
of the gravitational wave frequency, of the distance from the Sun, of the modulus of the declination and 
of the age such that respectively 50% and 90% of the detected sources have values lower than those. Only 
models giving N^et 1 ■ A are shown. 


Model 

Ndet 

A 

ho, 

ho[Hz] 

rso, 

7*90 [fcpc] 

\S\s,o, 

\0\so\deg\ 

ho, 

ho [Myr] 

1 

13.2 

288, 

473 

0.25, 

0.45 

18.3, 

49.8 

0.33, 

1.15 

2 

0.84 

490, 

954 

0.14, 

0.25 

18.3, 

54.4 

0.43, 

2.00 

4 

3.4 

190, 

198 

0.13, 

0.23 

20 .6, 

53.8 

0.25, 

1.48 

7 

8.3 

255, 

425 

0.19, 

0.41 

18.9, 

52.7 

0.34, 

1.43 

8 

0.82 

452, 

698 

0.13, 

0.26 

22.9, 

61.9 

0.53, 

2.05 




Figure 5. Model 1: distribution of the gravitational wave frequency and distance from the Sun for detected sources with the 
Virgo detector. 


first and second generation detectors. We use the sensitivity curves of Virgo detector and of 
a possible advanced Virgo detector, shown in Figs. m- 


5.1 Results for first generation interferometers 

Results are summarized in Tab. 0- The number of detected sources is parametrized by the 
fraction A ^ 1 of the total neutron star population made of GWDNS. Only the models giving 
N det > 1 ■ A are shown. For Model 1, for instance, we see that assuming A = 0.1, i.e. that 10% 
of all neutron stars are GWDNS, we would expect to have ~ 1 detection. More precisely, if 
Model 1 holds, the minimum fraction A m ; n we need to have at least one detectable source is 
7.6%. On the other hand, for Model 2 we would expect ~ 1 detection if all the neutron star 
population were made of GWDNS, i.e. A = 1, which could happen only if there were some 
mechanism suppressing the external magnetic field of a young neutron star, see Section 
ED Typically, detected sources rotate at high rate and are very near ( f gw > 200 Hz and 
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Figure 6. Model 2: distribution of the gravitational wave frequency and distance from the Sun for detected sources with the 
Virgo detector. 


r < 500 pc for most sources). Most detectable sources, moreover, are very young, with age less 
than ~ 2 Myr and then, some of them conld be detectable by ROSAT through their thermal 


emission, see e.g. 


Popov (2004). In such a short time interval the majority of GWDNS have 


moved for a rather short distance from their initial position, few hundreds of Megaparsec at 
the most, and then many detected objects were born in the solar neighbourhood and still 
belong to the local disk or to the Gould Belt. This explains why for all models about 50% 
of the detected sources have declination, in modulus, < 18°, which is the inclination of the 


Gould Belt, see Section rm If we consider only the loudest sources, however, the spread in 
declination is larger: half of the 1% loudest sources have a declination, in modulus, larger 
than ~ 26°. In Fig.© the gravitational wave frequency and distance distributions of detected 
GWDNS are plotted for Model 1; in Fig.© the same quantities are plotted for Model 2. 
For a given period distribution, the smaller is the mean cllipticity and the larger are the 
typical rotation frequencies of detected sources and the smaller are the distances. For Model 
1, for instance, the minimum gravitational wave frequency is ~ 80 Hz while it is larger than 
200 Hz for Model 2. Similarly, for Model 1 we have detected sources up to a distance from 
the Sun of about 8 kpc , even if, as already said, almost all are much nearer; for Model 2 the 
maximum distance is less than 0.7 kpc. No detection is obtained for Model 3. Fig. m refers 
to Model 4: note that the frequencies accumulate near 200 Hz which is the double of the 
maximum rotation frequency for that model; as the frequencies are on average lower than 
in Model 1, even if the mean cllipticity is the same, detected sources are on average nearer 
to the Sun. The minimum fraction we need in order to have at least one detected source is 
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Figure 7. Model 4: Distribution of the gravitational wave frequency and distance from the Sun for detected sources with the 
Virgo detector. 



Figure 8. Model 7: distribution of the gravitational wave frequency and distance from the Sun for detected sources with the 
Virgo detector. 


\ min = 29.4%. Models 5 and 6 do not give any detection. Perspective of detection for this 
kind of initial period distribution, would be obvioulsy even worse if the cut-off initial period 
were larger than 10 ms. Fig.® refers to Model 7. The frequency distribution is similar to 
that of Model 1, but a little narrower and no source is detected at frequencies larger than 
~ 1 kHz. Also the distance distribution is similar, but narrower, than that of Model 1. The 
minimum fraction A mm we would need to have -at least one detectable source is 12%. Fig.® 
refers to Model 8: from the figure and Tab.® we see that, as expected, detected sources are 
very near and emits signals at very high frequency. Like for Model 2, we could expect, at the 
most, one detection if all neutron stars were GWDNS. Model 9 produces no detection. In 
Tab.® we have reported, for each initial period distribution function and for three different 
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Figure 9. Model 8: distribution of the gravitational wave frequency and distance from the Sun for detected sources with the 
Virgo detector. 


Table 3. For each kind of initial period distribution and for three different 
values of the fraction A, in columns 2 — 4 the minimum allowed value of the 
mean ellipticity e in order to have at least one detectable source is reported; 
in column 5 there is the minimum allowed fraction if e = 10 —6 . The target 
Virgo sensitivity is assumed. 


pRo) 

A = 0.5 

A = 0.2 

A = 0.1 

Amin 

standard 

O 

1 

-4 

2.5 • lO" 7 

r- 

1 

O 

00 

0.076 

r- modes 

6.4 ■ 10“ 7 

- 

- 

0.27 

r-modes+fall-back 

2.4 ■ 10“ 7 

6.0 ■ 10" 7 

- 

0.12 


values of the fraction A, the minimum value of the neutron star mean ellipticity needed in 
order to have at least one detectable source; the minimum needed fraction A min , assuming 
e = 10~ 6 , is also shown. Only for A = 0.5 all the considered initial period distributions give 
at least one detection, with mean ellipticity in the range rv 2. -7- 6-10 7 . For A = 0.1 only 
the ’’standard” period distribution gives at least one detectable source (with e ~ 5 • 10 - '). 
Assuming e = 10” 6 , the minimum needed fraction varies between 7.6% to 27%. 

5.2 Results for second generation interferometers 

Results are summarized in Tab. ©• The detection rate increases by about two orders of 
magnitude respect to the case of detectors of the first generation. Fig. m refers to Model 
1. From it and from Tab. © we see that, while the frequency distribution is not heavily 
different from that concerning Model 1 for the Virgo detector, strong differences are evident 
for the other parameters. In particular, the distance distribution shows a large peak at small 
distances but also a second peak, with amplitude less than 50% lower, around 8.5 kpc, 
corresponding to the galactic centre. In any case, from the values of r 50 and ?’go it is clear 
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Table 4. Results for an advanced Virgo detector. The first column indicates the model, according to the 
classification of Tab. Ill: /50? ^50, |0|5O> ^50 and /90, ^90, |0|9o, £90 give the values of the gravitational wave 
frequency, of the distance from the Sun, of the modules of the declination and of the age such that respectively 
50% and 90% of the detected sources have values lower than those. The parameter A is the fraction of the 
total neutron star population made of GWDNS. Only models giving N^et > 1 • A are shown. 


Model 

Ndet 

X 

/so, fao[Hz] 

7*50 , rgo [kpc] 

050 , 090 [deg] 

*50, tg 0 [Myr] 

1 

1174 

295, 470 

3.8, 8.7 

1.03, 34.4 

0.17, 5.70 

2 

169 

378, 638 

0.51, 1.7 

19.7, 60.2 

3.65, 54.3 

3 

9.4 

518, 1245 

0.14, 0.36 

25.2, 66.2 

17.0, 314 

4 

240 

149, 195 

0.56, 1.35 

20.1, 60.2 

2.8, 33.4 

5 

19.5 

196, 199 

0.2, 0.44 

22.6, 59.0 

2.4, 189 

7 

635 

268, 449 

2.26, 8.5 

2.3, 43.8 

.38, 9.2 

8 

82.2 

334, 525 

0.45, 1.4 

22.6, 58.3 

4.3, 80 

9 

3.6 

440, 812 

0.14, 0.33 

25.4, 63.9 

51.3, 399 


0.09 

0.08 

0.07 

0.06 

>- 0.05 

D 

cr 

£ 0.04 
0.03 

0.02 

0.01 

0 

( 

Figure 10. Model 1: distribution of the gravitational wave frequency and distance from the Sun for detected sources with the 
advanced Virgo detector. 
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Figure 11. Model 2: distribution of the gravitational wave frequency and distance from the Sun for detected sources with the 
advanced Virgo detector. 
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Figure 12. Model 3: distribution of the gravitational wave frequency and distance from the Sun for detected sources with the 
advanced Virgo detector. 


that most of detected sources belong, in this case, to the galactic disk, and the contribution 
of the Gould Belt is small. This conclusion is reinforced by the value of 101 50 , which is around 
1°. However, if the loudest sources are taken into account, half of the 1% strongest emitters 
have declination larger, in modulus, than ~ 16.7°, comparable with the values we obtain 
for the Virgo detector, and belong to the local neighbourhood. Another difference, respect 
to the case of the Virgo detector, is that now we have a non negligible fraction of detected 
GWDNS with ages larger than ~ 5 Myr: this is an obvious consequence of the fact that 
due to the higher sensitivity we can detect sources emitting at lower frequencies, and then 
older. Fig. (HU) and Fig. m refer, respectively, to Model 2 and Model 3. As expected, the 
characteristic frequencies of the detected sources increase and the typical distances decrease, 
in order to compensate the smaller ellipticity. While for Model 2 we have still ~ 10% of 
the detected sources at r > 2 kpc, for Model 3 almost all the sources are in the solar 
neighbourhood. Another important difference respect to Model 1, and especially respect 
to the corresponding models for detectors of the first generation, is the characteristic age 
of detected objects which is now of the order of 100 Myr or even more. Fig. m and 
Fig. (HU refer respectively to Model 4 and Model 5. As already seen for the Virgo detector, 
for this kind of initial period distribution the detected sources have frequencies that tend to 
accumulate near the maximum possible value (200 Hz). This is particularly evident in Model 
5 for which the detection rate would rapidly decrease if the maximum possible frequency 
were even a little less than 200 Hz. No detection is expected for Model 6. Figs. mm 
refer, respectivley, to Models 7,8,9. For Model 7, as for Model 1, the detection distance 
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Figure 13. Model 4: distribution of the gravitational wave frequency and distance from the Sun for detected sources with the 
advanced Virgo detector. 




Figure 14. Model 5: distribution of the gravitational wave frequency and distance from the Sun for detected sources with the 
advanced Virgo detector. 


arrives beyond the galactic centre and this produces the secondary peak in Fig. (jT5l) . Most 
of the detectable sources are near the galactic plane but, as for Model 1, the loudest sources 
are more spread: half of the 1% strongest emitters has a declination, in modulus, larger than 
16.6°. The distance distributions of Models 7-9 are rather similar to those we find for Models 
1-3 while concerning the frequency distribution the main difference is that the maximum 
frequency is always less than 1000 Hz, due to the minimum initial period equal to 2 ms. In 
Tab.© we have reported, for each initial period distribution function and for three different 
values of the fraction A, the minimum value of the neutron star mean ellipticity needed in 
order to have at least one detectable source; the minimum allowed fraction A min , assuming 
e = 10~ 6 , is also shown. A value A = 10% would allow us to detect sources with mean 
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Figure 15. Model 7: distribution of the gravitational wave frequency and distance from the Sun for detected sources with the 
advanced Virgo detector. 




Figure 16. Model 8: distribution of the gravitational wave frequency and distance from the Sun for detected sources with the 
advanced Virgo detector. 


ellipticity as small as ~ ICE 8 . Even for A = ICE 3 we could still have a detection, if the mean 
cllipticity is e = 8.6 • 10~ 7 and if the period distribution were described bt the ’’standard” 
model. If I — ICE 6 , the minimum needed fraction varies between 8.5 • ICE 4 and 4.2 ■ ICE 3 , 
depending on the initial period distribution. 


Table 5. For each kind of initial period distribution and for three different 
values of the fraction A, in columns 2 — 4 the minimum allowed value of the 
mean ellipticity e in order to have at least one detectable source is reported; 
in column 5 there is the minimum needed fraction if e = 10 -6 . The advanced 
Virgo sensitivity shown in Figs. 03 is assumed to hold. 


p(Po) 

A = 0.1 

A = 0.01 

A = 0.001 

min 

standard 

10 -8 

6 ■10-® 

8.6 • 10" 7 
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00 

r-modes 
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Figure IT. Model 9: distribution of the gravitational wave frequency and distance from the Sun for detected sources with the 
advanced Virgo detector. 


6 POSSIBLE MECHANISMS ALLOWING THE EXISTENCE OF A 
POPULATION OF GWDNS 


GWDNS could be the result of an evolutionary path different from that of pulsars or conlcl 
belong to the same population of pulsars, but with a very low dipolar magnetic field. Such a 
low magnetic field could be produced by several mechanisms some of which are now shortly 
discussed. 


According to 


Geppert et al 


(1999), post core-collapse hypercritical accretion can sub¬ 


merge any initial magnetic field of the newborn neutron star, because the ram pressure at 
the star surface is much larger than the magnetic field pressure. Later diffusion of the mag¬ 
netic field could produce a delayed switch-on of a pulsar, with timescale up to millions of 
years. There are, however, many possible factors which could reduce or even suppress the 
fail-back, like a high kick velocity or the injection of rotational energy in the highly dirty 
environment of the new-born neutron star. It is not possible at present to estimate which 
fraction of newborn neutron stars undergoes a phase of hypercritical accretion. As most 
detected objects have, in our simulations, ages below 10 Myr (particularly in the case of 
the inital Virgo detector), this mechanism could be efficient in producing a population of 


GWDNS. 


If the magnetic field of neutron stars is associated with the alignement of the mag¬ 
netic moment of neutrons in the stellar crust, the orientation of the magnetic field will 
be influenced by the coupling between intrinsic spin and gravitational field of the star. 
This coupling induces the alignement of the magnetic axis respect to the rotation axis 
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(Zh ang. Cheng & Wulll998jh In such a situation, IZhang. Cheng & WiJ (jl998li have shown 
that the timescale of alignement is inversely porpoprtional to the magnetic held strength 
and proportional to the pulsar initial period. As a consequence, the inclination angles of 
weak field-short initial period pulsars should concentrate toward lower values. The relation 
they find for the time evolution of the inclination angle a is 

/ Bin VV B 


F[a(t)] = F[a(0)] - 0.023 


V10 12 G) VlO 12 G 




(26) 


where F[a\ = a — | sin 2a and B in is the star inner held. 

A third possible mechanism is the decay of the neutron star magnetic held. As it is well 
known, the time evolution of the magnetic held in isolated neutron stars is still a controversial 
issue. A hrm point is that the observation of radio pulsars rules out decay times less than 
~ 10 Myr but larger decay times are still possible. The faster is the decay of the magnetic 
held and the larger is the fraction of GWDNS in the neutron star population. 


7 CONCLUSIONS 

In this paper we have described the Monte Carlo simulation of a population of asymmetric 
neutron star which evolve through the emission of gravitational waves and have estimated 
the detectability of the emitted continuous gravitational signals. As far as we know, a quan¬ 
titative study of this kind of sources has not been done before. Some important features of 
this work are: neutron star cllipticity is not considered a constant, but distributed according 
to an exponential law; the presence of the Gould Belt is taken into account; the number 
of detectable sources is estimated assuming a realistic data analysis method is applied (in 
fact, the analysis method we will apply to the data of the Virgo detector, once they will 
be available). We have performed the simulation for values of the neutron star mean ellip- 
ticity between 10 -6 and 10 -8 and for three different models of the initial rotation period 
distribution corresponding, respectively, to a ’’standard” pulsar population, to a population 
of neutron stars spun-down by r-modes and to a population subject, after the supernova 
explosion, to the combined effect of r-modes and matter fail-back. The number of detections 
is estimated using the Virgo target sensitivity and a possible sensitivity curve of an advanced 
Virgo detector. No large difference is expected for LIGO interferometers because very few 
sources are likely to be detected in the low frequency range, say below 40 Hz, where the tar¬ 
get sensitivity of Virgo is much better than that of LIGO. We do not have any clue of what 
the number of GWDNS could be, then the number of detectable sources is parametrized by 
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the fraction A of the total neutron star population made of this kind of objects. For detectors 
of the first generation, we need a fraction A ~ 8% to have one detectable source, for e ~ 10 -6 . 
Smaller values, e ~ 2 -t-5 • 1CT 7 , require a fraction A as high as ~ 50%. The expected number 
of detections for advanced interferometers is about two orders of magnitude larger than that 
we find for initial detectors. This means that, assuming a value for A, for advanced detectors 
the minimum detectable mean cllipticity is much smaller or, viceversa, for a given value of 
the mean cllipticity, the minimum needed fraction of neutron stars being GWDNS, in order 
to have at least one detectable source, is much lower. We find that a value of A as small as 
0.001 is enough to have one detectable source for high mean cllipticities, while if A ~ 10%, 
the minimum needed mean cllipticity is of the order of 1CD 8 . 

In principle, GWDNS as old as the Galaxy could contribute to the detected popula¬ 
tion, because for many objects the slow gravitational spin-down does not bring the signal 
frequency below the detector low frequency cut-off even for age of the order of 10 10 yr. 
However, in practice detectable sources are often very young, especially in the case of inter¬ 
ferometers of the first generation, with ages of few Megayears or less. For many models, most 
detectable sources are very near to the Sun, a few hundreds parsec at most, and the Gould 
Belt contributes significantly to the detectable population. If these sources are really located 
in the solar neighbourhood, some of them could be detected in the electromagnetic band 
through their thermal emission. The distance ”reached” by interferometers of the second 
generation is much larger than that of initial detectors: even sources located in the galactic 
centre can be detected, for some models, but strongest sources are expected to be near and 
then distributed rather homogeneously around the Sun. This, and the poor knowledge we 
have about typical frequencies and spin-down rates of GWDNS, are the reason why a blind 
search is needed in any case. As it is by now well known, a blind search can be performed 
only using non optimal data analysis methods with reduced sensitivity respect to optimal 
ones (typically a factor of 2-3) but which require a reasonable, although large, amount of 
computing resources. More precisely, the larger is the available computing power and the 
wider can be the search, covering a larger portion of the source parameter space. 

This work has shown that GWDNS are potentially interesting sources of gravitational 
radiation. Clearly, the very poor knowledge we have about parameter values, in particular 
the ellipticity, prevents us from giving firm estimations of the expected number of detections. 
By the way, even for interferometers of the first generation, already taking data or being on 
the point to start data collection, it is possible, although rather unlikely, that some detection 
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takes place. Possibilities of detection are obviously much better for interferometers of the 
second generation, which could be able to catch several sources. Indeed, gravitational wave 
astronomy would be the best tool to study the properties and demography of GWDNS. 


APPENDIX A: OTHER MECHANISMS COMPETING WITH THE 
GRAVITATIONAL EMISSION OF GWDNS 

In the previous sections we have considered GWDNS as isolated asymmetric neutron stars 
with a magnetic held strength low enough so that the gravitational spin-down is dominant 
respect to the standard dipolar magnetic held spin-down. However, as already anticipated 
in Section 2, a neutron star can pass through different evolutionary stages in which the spin- 
down is due to other mechanisms, like those connected to the interaction of the star with 
the interstellar medium (ISM). In this section we shortly remind the standard classification 
of isolated neutron stars evolutionary stages and try to understand if mechanisms different 
from the dipolar held spin-down can overcome the grav itational spin-down in GWDNS. 


According to the standard classification (Popov et al. 


2000j), basically three different 


evolutionary stages can be outlined for neutron stars. In the ejector phase the gravitational 
energy density at the accretion radius is lower than the energy density of the momentum 
outhow produced by the rotating magnetic held; neglecting the possible emission of grav¬ 
itational radiation, the spin-down is due to the rotating dipolar magnetic held. A neutron 
star in this phase can appear, but not necessarily, as a pulsar. This phase, depending on 
the strength of the held and on the star velocity, can last even for the whole neutron star 
life. The propeller phase starts when the gravitational energy density at the accretion radius 
is larger than the energy density due to the momentum outhow produced by the rotating 
magnetic held. This condition is verihed when the rotation period is larger than 
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is the accretion radius. The accretion rate on a neutron star, in the Bondi theory, is 
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where n is the particle density in the ISM (~ 1 cm 3 ) and v is the neutron star velocity. 
Another characteristic length which we will use in the following is the Alfven radius 
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In the propeller phase the neutron star can start to interact with the interstellar medium: 
matter can fall toward the neutron star up to the Alfven radius where the magnetic energy 
density balances the kinetic energy of the falling matter. The co-rotating magnetosphere will 
prevent the accretion of matter on the stellar surface, unless the gravitational acceleration at 
the Alfven radius is larger than the centrifugal acceleration and the Alfven radius is smaller 
then the accretion radius. The latter condition implies that the neutron star velocity satisifes 
the inequality 
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The former condition translates in a second critical period 
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For periods P > P A we are in the accretor phase where the matter can fall up to the stellar 
surface. Both in the propeller and in the accretor phase the neutron star spin-down is not 
given by Eq. ©, which holds in the ejector phase, then the condition to have a GWDNS, 
given by Eq.(|%|). must change. From Eqs. muD we have that the condition the magnetic 
field must satisfy in order to bring the neutron star in the propeller phase is 
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10 km,s~ l , 

Here and in the following equations the neutron star is supposed to have a radius R = 10 km. 
Assuming the neutron star is in the propeller phase, let us see which is the further condition 
we need on the magnetic field so that the spin-down is still dominated by the emission of 
gravitational waves. In the proyell er phase the neutron star looses its rotational energy to 
the ISM and spins-down at a rate (Popov et al. 2000) 
fMr 2 c 
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is the corotation radius , 
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u k {r A ) = 


GM 


(A10) 


is the local keplerian velocity at r A and £ is a numerical factor depending on the geometry of 
the accretion flow, which for spherical accretion, typical of isolated neutron stars, is ~ 10“ 3 . 
Using Eq. (ED and Eq. (1A8I) and imposing that Y' = ■^ S2S - > 1 we find: 
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For instance, assuming e = 10 , v = 10 fcms , / = 25 Hz we have B crit:Prop ~ 10' G 
and B' critprop ~ 6 • 10 4 G, that is, even if the magnetic field is so low that the neutron 
star moves to the propeller phase, the gravitational spin-down is still dominating, unless the 
magnetic field is nearly vanishing. This conclusion is even strengthened if we consider that 
most neutron stars have typical velocity much larger than 10 kms~ l and that the accretion 
rate given by Eq. (HI is only an upper limit and is probably much lower in real situations. 
From Eas. (IA6IA3ll we obtain also the condition on the magnetic field in order to enter in 
the accretor phase: 
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If this condition is satisfied, Eq. HS l which must be also satisfied to be in the accretor 
regime, gives v < 6 • 10 3 kms _1 , and then is also verified. Now we compare the spin-down 


rate of a neutron star in the accretor stage with that due to grav itational radia 


ion. The 


spin-down rate in the accretor stage can be written as (Konenkhov & Popov 1997 1 
B 2 R 6 


- firr kf 


I<i 


turb 


(A 13) 


where the first term in the righ-hand side describes the braking moment of the forces due to 
the possible turbulization of the ISM, k t is a dimensionless costant with value 1 and 
is a random term which superimposes to the first one and can produce, on small time scales, 
both spin-down and spin-up of the neutron star. Imposing that Y" = > 1 we End our 

last condition on the magnetic field strength in order to have a neutron star dominated by 
the emission of gravitational waves: 
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If Ea. (IA12|) is verified also Ea. (IA14j) is satisfied, at least if we limit ourselves to consider the 
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frequency range where detectable (in principle) sources can be placed, that is / > 25 Hz. 
As a conclusion, we can say that, unless the magnetic held is nearly vanishing, a neutron 
star which is not in the ejector phase will be dominated by the emission of gravitational 
radiation. 
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